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Giant oscillations of the density of electronic states and the differential conductance of a 
superconductor-ferromagnet-superconductor structure are predicted for the case that the exchange 
$_l ' energy of the interaction between the electron spin and the spontaneous moment of the ferromagnet, 

7o, is less than the superconductor energy gap, A (Io < A ). The effect is due to an extremely large 
degeneration of the energy level e — Io (e is the electron energy measured from the Fermi-energy) 
if the superconductor phase difference tp is close to odd numbers of tt. These quantum interference 
effect persists even in long ferromagnetic bridges the length of which much exceeds the "magnetic 
length" Tivf/Iq for the ballistic case and ^/ hD / Io for the diffusive one (D is the electron diffusion 
constant). The predicted effect allows a direct spectroscopy of Andreev levels in the ferromagnet as 
well as a direct measurement of the exchange energy,7o, of the interaction of the electron spin with 
' j^j ■ the spontaneous moment of the ferromagnet. 
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INTRODUCTION 



In recent years much attention has been paid to the conductance of mesoscopic superconductor-normal conductor- 
superconductor systems (S/N/S heterostructures) which is sensitive to the phase difference between the supercon- 
ductors (see, e.g., the review paper by Lambert and Rimondi and references there). This effect takes place both 
in ballistic and diffusive samples due to a quantum interference caused by Andreev reflections of quasi-particles at 
two (or more) N/S interfaces which imposes the phase of the superconducting condensate on the quasiparticle wave 
q , function in the normal conductor. 

' In normal conductor-superconductor heterostructures, electronic elementary excitations which freely propagate 
l— ~ 1 , inside the non-superconducting metal can not penetrate into the superconductor if their energy e (measured from the 
I ' Fermi energy) is less than the superconductor energy gap A. A correlated transferring of two electrons accompanied 
, by their pairing inside the superconductor is the only mechanism that provides a direct transmission of the charge 
into the superconducting condensate that is the ground state of the superconductor. This 2-electron transmission 
may be considered as a scattering process under which electronic excitations inside the normal conductor undergo 
an electron-hole transformation at the boundary with the superconductor. This scattering (known as the Andreev 
reflection) couples the incident electron (hole) and the reflected hole (electron) in such a way that their spins are 
oriented in opposite directions and their energies (±e) are symmetrically positioned with respect to the Fermi energy 
£f ("Andreev hybrid"), as shown in Fig^ As a result of the Andreev reflection the electronic excitations reflected 
at the normal conductor - superconductor interface pick up the phase <f> of the superconductor order parameter 
A = |A|expi0 and keep memory of it inside the normal conductor. Such a 2-electron superconducting correlation 
G ' persists inside the normal conductor at the distance calL e from the superconductor, 

C e = min (h/Ap, JhD/Apv^) (1) 

o: • ; 

where Ap = p e — Ph ~ s/vf, vp is the Fermi velocity and D is the diffusion coefficient. At larger distances from 
the superconductor the destruction of the phase coherence arises due to the difference between the momenta of the 
electron (p e = pp + s/vf) and hole (ph — Pf — s/vf) components in the Andreev hybrid (the typical value of the 
energy is s ~ kT, where T is the temperature and k is the Boltzmann constant). The other peculiar feature of 
the Andreev hybrid is that the electron and hole spins have opposite directions. This spin flip does not change the 
interference pattern of the non-magnetic metal because all energy levels are twice degenerated with respect to the 
spin direction. In ferromagnets, however, this degeneracy is lifted due to the interaction of the electron spin with the 
spontaneous moment of the ferromagnet (below we refer to it as the exchange- interaction energy To), and electrons 
with opposite directions of the spins occupy different energy bands (see Fig|2): 

,2 



£<r(p) = p72m + oT (2) 



where p is the electron momentum, m is the electron mass, a = ±1 for the electron spin up and down. In a 
ferromagnetic metal-superconductor heterostructure, the change of the spin direction of the incident electron (hole) 
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FIG. 1: Schematic representation of the Andreev reflection: the incident electron with energy e (measured from the Fermi 
energy ef) and the spin up is reflected back as a hole with the energy — e and the spin down; the incident electron and the 
reflected hole momenta are p e = Wp 
respectively. 



2me and p e 



2me (pf is the Fermi momentum, m is the electron mass) 




FIG. 2: Energy bands for electrons with opposite spins 



under the Andreev reflection shifts the reflected hole (electron) into the other energy band that causes an additional 
difference dp = Ia/vp. The latter drastically decreases the penetration length L by orders of the magnitude |6^| 
that influences the quan tum interference affecting many propert ies of ferromagnet /superconductor heterostructures 
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a shortening of the proximity effect has been actually observed in ma gnet ic materials Refs. jo, |3, llii LLJ, On 
the other hand, measurements carried out in recent works [32t 133 . 134L 35, 36] demonstrate a long-range proximity 
effect in magnetic materials that is in an obvious contradiction with the above general considerations. As shown in 
papers |37l l38| such a long-range proximity effect arises in F/S heterostructures in which the ferromagnetic part is 
magnetically inhomogeneous. In this case a triplet su per conducting correlation persists on a length typical for non- 
magnetic materials j22i Ha, Uls EH C3 El EE EE E3 ■ A long-range proximity effect of another type occurs 
in F/S heterostructures in which the ferromagnetic part is magnetically homogeneous but the exchange-interaction 
energy I is less than the superconductor energy gap A. In this case the difference between the electron and hole 
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momenta of Andreev hybrids with energy e = Jo do not depend on Jo (see the section^ and hence their penetration 
length is equal the one for non-magnetic conductors C e |48|, |4!j It results in giant oscillations of the conductance of an 
S/F/S Andreev interferometer with a change of the phase difference between the superconductors when the voltage 
V applied to the ferromagnet is close to 2I /e |4^; for the case of F/S structures the subgap conductance exhibits a 
peak below the superconducting energy gap A at eV = Iq < A 0| 

It is important for applications that the both above-mentioned types of the long-range proximity effects allow to 
create superconducting quantum interference devices with a ferromagnetic normal metal junction of an anomalous 
large length (a device of such a type was suggested in [5tT]1. 

In this paper we consider the density of states, the differential conductance and the current-voltage characteristic 
(CVC) of an superconductor-ferromagnet-superconductor heterostructure (S/F/S) of the Andreev interferometer type 
in which the ferromagnetic part is separated from the reservoir of normal electrons with a potential barrier of low 
transparency |63, t r <l, (see Fig|3J|. For the case that the exchange- interaction energy Iq is less than the supercon- 
ductor energy gap A we predict a sharp peak in the density of states at energies close to Iq and sharp dependence 
of the differential conductance G(V, ip) on the applied voltage V that is accompanied with giant oscillations of the 
density of states and the conductance with a change of the phase difference ip between the superconductors. For 
the geometry under consideration the differential conductance is proportional to the density of Andreev states that 
allows their direct spectroscopy with electrical measurements. On the other hand, the predicted effect allows to find 
the electron spin - ferromagnet moment exchange interaction energy Iq with a direct electric measurement because a 
sharp and high peak in the differential conductance takes place exactly at V = 2Iq. 

In section B we qualitatively explain the phenomenon under consideration. In sections B and B we present analytical 
and numerical calculations for the density of states and for the current and the differential conductance in ballistic 
S/F/S structures, respectively. In section „ we calculate the density of states in a disordered S/F/S structure using 
the Gutzwiller path-integral approach. 



In this section we qualitatively show that a change of the superconductor phase difference ip can result in gi- 
ant oscillations of the conductance of a superconductor-ferromagnetic conductor - superconductor structure if the 
exchange-interaction energy Iq of the ferromagnet is less than energy gap A of the superconductor. The effect arises 
due to a resonant transmission of quasi-particles from the normal reservoir through Andreev states in the ferromagnet 
which macroscopically concentrate near the exchange-interaction energy Iq if the superconductor phase difference ip 
is close to an odd number of tt and the voltage V applied to structure (see FigOJ) is close to 2Iq/c (Iq < \A\). A 
qualitative explanation of the effect is as follows. 

Under Andreev reflection at an F/S interface the spin directions of the incident and reflected quasi-particles are 
opposite and hence the longitudinal momenta (parallel to the cc-axis which is perpendicular to the F/S interfaces) of 
the electron and the hole are, respectively, (de Jong and Beenakker 52] ) 



where pp is the Fermi momentum, pj_(n) = (0, Tin y /d v , hn z /d z ) is the quantized transverse quasi-particle momentum 
parallel to the F/S interfaces, (d v and d z are the transverse sizes of the ferromagnetic section, n = (0, n y , n z ); n yi n z = 
0,1,2,..., |n| <N±_ are the transverse mode quantum numbers) assuming a hard wall confining potential with the 
number of transverse modes N± w S/Xp inside it, S — d y d z is the cross section area of the ferromagnet (equal to the 
area of the F/S interfaces), = 2-kTi/pf is the Fermi wave length. From Eq. 10 it follows that in contrast to the 
non-magnetic case, near the Fermi level (e « 0) the electron and the hole momenta in the ferromagnet are different, 
and for a large enough I (usually Iq is greater than the Thouless energy ) the interference effects are absent due to 
the destructive interference. This fact demonstrates the conflict between superconductivity and magnetic ordering in 
S/F/S structures. 

However, interference effects in the ferromagnet can exist albeit at some finite voltage V applied between the 
reservoir and the superconductor. Our argument starts from a description of the electron transport in terms of 
resonant tunneling through quantized energy levels of the ferromagnet mesoscopic part of the system shown in Fig. 
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FIG. 3: Superconductor-ferromagnet-superconductor heterostructure of the Andreev interferometer type. Thick lines indicate 
potential barriers of low transparency, t r <C 1, separating the ferromagnet from the reservoir of normal electrons; voltage V is 
applied between the reservoir and the superconductor. 



Taking into account the amplitude of the Andreev reflection at a normal conductor - superconductor interface |5- 



one easily finds the semiclassical quantization condition for an S/F/S system in the absence of potential barriers at 



where pi^o-(e) and pn,o-(e) ar e the longitudinal momenta determined by Eq. (J3J; L is the distance between the F/S 
interfaces; ip = <pi — <f>2 is the phase difference between the superconductors 1 and 2; I — 0, ±1, ±2, .. is the longitudinal 
quantum number; while writing equations Eq.© and Eq.© we assumed e <C |A| for simplicity's sake. Dispersion 
equation Eq.© determines Andreev levels £; !n which can be shifted with a change of ip by means, for instance, of 
an external magnetic field. When an Andreev level is lined with the energy of electrons injected from the reservoir 
of normal electrons, the resonant transmission of electrons through the ferromagnetic section takes place that causes 
an increase of the system conductance. Typically each Andreev level is only 2-fold (spin up and down) degenerated, 
and hence simultaneously only 2 electrons can be resonantly transmitted through this energy level. It means that the 
amplitude of oscillations of the differential conductance G(V) = dJ/dV (J is the dissipative current current, see Fig. 
© with a change of the superconductor phase difference ip is of order of e 2 /h. The situation drastically changes when 
the superconductor phase difference is equal to an odd number of ir. In this case the energy level e = Iq is highly 
degenerated since all N± transverse modes are simultaneously at this level (if ip = 7t(2Z + 1) the dispersion equation 
Eq.© is satisfied with e — Iq for all transverse modes because p^ = p^ 1 ' for any n, see Eq.© and Eq.©). Therefore 
the resonant transmission occurs simultaneously in all transverse modes when the superconductor phase difference 
is equal to odd numbers of it and the applied voltage takes such a value that eV/2 is equal to the exchange energy 
Jo, thus producing a giant conductance peak. The width of the peak in the dependence of the conductance on ip at 
V = 2Iq / e is of the order of the single-electron transparency t r of the potential barriers between the ferromagnet and 
the electron reservoir, and the width of the peak in the ^-dependence of the conductance at <p = 7r(2Z + 1) is ~ t r Eo 
(Eq = Hvf/L is the distance between neighboring energy levels); the amplitude of the peak is of the order of N±e 2 /h 
reflecting the above-mentioned N i -fo ld degeneracy of the resonant level. When this result is compared with that for 
a non-magnetic normal part (see |55l l56j 1 it is apparent that the exchange interaction of the electron spin with the 
ferromagnet spontaneous momenta destroys the giant oscillations of the conductance on the Fermi level (that is at V 
close to zero) shifting them to the energy range eV/2 « Io where the giant conductance oscillations are restored even 
if/ >£o. 

Giant conductance oscillations for a diffusive ferromagnet in the case that the temperature T satisfies the inequality 
t T E^? <C kT <C = TiD/L 2 were considered in paper |4g|. In section, we show that the high degeneracy of the 



ta = exp (i7r/2 + (f>) 



(4) 



F/S interfaces to be 
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energy level e = Iq at ip = ir is not lifted in diffusive S/F/S structures that results in a sharp peak in the density of 
states and hence in giant conductance oscillations 



DENSITY OF STATES IN BALLISTIC S/F/S STRUCTURES WITH I < |A| 

In this Section we find the dispersion equation and the density of Anreev states in a ferromagnetic conductor- 
superconductor- ferromagnetic conductor structure in which potential barriers are present at the ferromagnet- 
superconductor interfaces. 

In order to find Andreev energy levels e a n l inside the ferromagnetic conductor we use the Stoner model with the 
exchange interaction energy Iq so that the Bogolubov-de Gcnnes equations are written as follows 52] : 



I [H + o-I(r) - ej u +a + Av. a = 
\ A*u +a - (H - al{r) + e)v- a =0 

where Hq = p 2 /2m— ep\ the superconducting energy gap A(r) and the ferromagnet exchange energy I(r) have non-zero 
values in complementary space regions, A = const ^ 0, I = in the superconductor and A = 0, I = Iq = const ^ 
in the ferromagnet. We assume Andreev reflections of quasi-particles at the F/S interfaces to be accompanied by 
normal reflections, that is the scattering process at an F/S interface is described by a 2 x 2 scattering matrix |58Ll59j 
as follows 68] : 

V ~ r A r N exp(-i?7) J 

where 

•\/|io| 4 + 4|r | 2 sin 2 ^ 
exp(-iip e ) - |r | 2 exp(zi/> e )' 

e*. - . If', ; (8) 

In Eq.ljJJ) the probability amplitudes of Andreev (ta) and normal (rjv) reflections are written as follows: 

i\to\ 2 exp(i(f>) 2ros'mip £ 
ta = =; r N - 



^Itol 4 + 4|r | 2 sin 2 Ve ' ^M 4 + 4|r | 2 sin 2 ^ e 



where cj> is the phase of the superconductor energy gap A = |A| exp(z</>); ro and to are the probability amplitudes for 
an incident electron to be reflected back and to be transmitted through the interface in case that the conductors on 
the both sides of the interface are in the normal state; this scattering arises due to a potential barrier at the interface, 
mismatch between the Fermi velocities of electrons of the conductors or of their effective masses, and so on. In 
general case this scattering is described by the scattering matrix 

,„ n ro] esp(W) i|fo| V \ rQ f + \ to f = 1 

\i\to\ | r \ exp(-i$) J 

According to Eq.© the electron- and hole-like components of the wave- function in the n-th transverse mode inside 
the ferromagnet is 



u a (x,y) = (<$**®* +b$e-*&*) sinfc ± (n)y 

v a (x,y) = (aWe- ik » )x + ftWe"^*) smk ± (n)y (9) 

Here and [a^ and bffl] are the probability amplitudes for free motion of electrons [holes] forward and 
backward, respectively, in channel n; kn' h ^ = Pn^/h (see Eq.@) is the electron (hole) longitudinal momentum for 
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an electron (hole) in the n-th transverse mode and its spin direction a (a = ±1), while x and (y, z) are longitudinal 
and transverse coordinates in the sample, respectively; here and below {n} stands for the summation over all the 
open modes |n| < Nj_. 

Matching the wave- functions Eq.© at two nonequivalent F/S boundaries with the use of the scattering matrix 
Eq. (J7J results in a spectral function of the form 



D n } ( e ) = cos l P~ -\ r N\ \ r N I cos f+ + \ r A ] I \ r< A I cos <P 



(10) 



Here </?-(e) = (fen — kn)L + x+ an d <P+(e) — (fen + kn)L + fx+ where L is the length of the ferromagnet, 
X+( £ ) = Xi + X2', A*+ = ? /i + ^i + T)2 + J?2, and <p — <f>i — 4>2 is the phase difference between the superconductors; the 
indices 1 and 2 indicate the F/S boundaries. The Andreev discrete energy levels ; of the system are determined by 

solutions of the equation (e) = 0, that is 



({krijj 



ki h l)L - 



X+ 



\r$ | |r£> | cos + fc£> )L + » + ) - \ \rf | cos <p 



(11) 



In the case that the ferromagnetic part of the structure is coupled to a reservoir of normal electrons through a 
potential barrier of a low transparency t r <C 1 (see Fig. |3J), the density of Andreev states inside the ferromagnetic 
conductor can be written as follows: 



(12) 



<7=-l {n} 



where V is the volume of the ferromagnetic part, is the density of Andreev states at a fixed transverse mode 
quantum number n and the spin projection a: 



7T ' ■ 



t r Eo 



1 



(e - e'Y + (t r E y 



(13) 



and 



(14) 



For convenience sake of further analytical calculations we write the density of Andreev states in the following 
form 



i 



dD<£\e) 



ds 



5(D£\e) 



(15) 



where (e) is defined by Eq. ijTU|) . 

To find the density function v^{e) we use the method developed in Ref. [6l]]. As dip-/de w ±l/E n (here = 
%v n /L and w n = yV^ — p^(n)/m) the factor dDn\s) / de is a trigonometrical function of ip±, as well as I?i' T ' > (see 
Ea.pOJl). and it is productive to expand i/„ into Fourier series in ip± and write it as follows: 



i— — OQ S — — QO 



A"£ exp(is< y £>_(e) + ifey> + (e)); 



(16) 



A* 



(27T) 



de 



5 f- ^ (¥>+>¥>-)) exp(-is<^_ -ikip + ) 



(17) 
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FIG. 4: Normalized density of electronic states inside the ferromagnetic conductor in a ballistic S/F/S structure for Iq = 0.5i?o 
r (i) _ 1 r {2 



and rty = 0.1, r$ = 0.05 and t r = 0.1; Vq = 8N ± /(nE );. Fig. A shows the dependence of the density of states on energy e 



at tp — n (full line), tp = 1.17T (dash line), and tp = 1.2n (dash-dashed line); Fig. B shows giant oscillations of the density of 
states with a change of the superconductor phase difference tp at e = Jo. 



As shown in Ref.|5a| the main contribution to the state density function Ea. (|16fl is of the terms A"q which, after 
integrating Eq. ((TBjl with respect to <p- , is written as 



.4 



8e ) (2tt) 2 



dip + {exp [isipi(<p+)] + exp [-istp 1 ((p + )]} 



(18) 



where 



V?i(v?+) = arccosdr^Hr^l cosv?+ - \r%> \ \r { ^ ] \ cos ip), < < tt 
Inserting Ea. (|18|l into Ea. (1161 (other terms in Ea. l|16|l with k ^ are negligibly small, see Ref . p3(| ) one find v° to be 



IT , , / , " A + 

"- (e) = l e: + ~dT 



2 E 7^2 rV{exp[i S (^(^) + ^-(^))]+exp[i S (-^(^)+^-(£))]} (19) 



Using Ea. (|19|l . Ea. (|13|l and Ea. (|12f> one finds the the density of electronic states in the ferromagnet to be equal to 



V de 



"dtp, 



2ir l-cos[2(e + aI )y/E + X +(£)+<r'ipi{<P+)} + 2t 2 r 



(20) 



Numerical results for the density of states based on Eq. (|20f> are shown in FigQ] The sharp peak in the dependence 
of the density of states on energy corresponds to the A^-fold degeneracy of the energy level e = Iq in the case that 
the superconductor phase difference ip is equal to odd numbers of 7r, as was explained in the previous section. 



NORMAL CURRENT AND CONDUCTANCE OF S/F/S STRUCTURES WITH J < |Aj 

In this Section we find the differential conductance and the current through the ferromagnet in the geometry of 
Fig. |21 for the case that the ballistic ferromagnet is weakly coupled to a normal electron reservoir (that is t r < 1) 
with a voltage drop V to be applied between the reservoir and the superconductor. We assume that there are 
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('12') ( 1 2) 

Andreev and normal reflections at F/S interfaces 1 and 2 with the probability amplitudes r A : and r y N ' , respectively 
(|r^| a + |r<^| a = l). 

According to the Landauer-Lambert formula |57j the current from the normal electron reservoir into the supercon- 
ductor is written as follows: 



J = ££(/ ( £ -f)-/ o(£ + f))^( £ ); 



(21) 



(n) 



where fo(s) is the Fermi distribution function, is the probability for an electron approaching the ferromagnet 

along the lead to be reflected back into the reservoir as a hole. 

As we are interested in the resonant transmission of quasi-particles through Andreev levels l inside the ferromag- 
net we use the resonant Breight-Wigner formula to write Ra{s) in the form 



(22) 



For convenience sake we represent Ea.()22|l as follows: 

+i 



[{e-e')/E n f+tl 



(23) 



<r=-l {„}' 

Inserting Eq. l|19|l into Ea. l|23[l one finds the electron- hole transmission coefficient -Ra(e) to be equal to 

l l 



2N ± x - x - 

r a (s) = —± J2 E 



<j=-l a'=-l J 1 y 



y de 



2fi 



, 2tt l-cos[2(e + oJ )tf/^+X+(e)+^i(¥>+)] + 2 *? 
Using Eq.J23J an d Ea. lj2*T|) one finds the differential conductance G(V) = dl/dV to be 

2e 2 N, ' 



G{E) 



nil 



12 E E 

<r=— 1 cr'=-l 



r dip + 



i 2T Jo 27r 



(24) 



2/ 



1 - cos [{eV + 2al )y + a'<pi(<p+)] + 2t 2 r 
where V — V/En and In = Iq/ E, and the current to be 



J = 



8eN ± E t r 



7Th y/T+Vr^h y'Jo 



E 



dy r dip, 



2tt 



(25) 



arctan 



V^TT^, { (eV/2 + I )y Vl {<p, 
tan — h o — 



En 



arctan 



tan 



-eV/2 + I )y gifo 
h cr — 



En 



(26) 
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FIG. 6: Normalized current- voltage characteristics for phase differences tp = n (full line), tp = I.l7r (dashed line), and tp — 1.2n 
(dash-dotted line) shown for r$ = 0.1, rff = 0.05, U = 0.1, and I = E ; J = 2eN±E /h 



Numerical results for the differential conductance and the current based on Ea. (|25|l and Ea. (|2fi[) are shown in Fig. 
03 and Fig. El They demonstrate a high sensitivity of the differential conductance and the non-linear current-voltage 
characteristics to both the superconductor phase difference tp and the applied voltage V. 

At low voltages, far from 2I /e, we have a resonant tunneling of quasi-particles through separate Andreev levels, 
and the conductance and the current are low. When eV/2 « Iq and tp = n(2l + 1) (/ = 0, ±1, ±2, ...) Andreev 
levels concentrate near Iq as can be readily seen from Ea. H10(l . and we have simultaneous resonant transmission of 
quasi-particles through the whole number of N±_ states resulting in a high peak in the conductance and a sharp jump 
in the current. When tp deviates from an odd number of tt, the number of Andreev levels concentrated near Iq is 
decreasing that results in a decrease of the sensitivity of the conductance and the current to the voltage. 

As is evident from Ea. (|ll|l . in case r^J ^ rjy Andreev levels are repelled from the level Iq. It results in a splitting 

(proportional to <5rjy = \r*jp — r$ | , see |5dj ) of the conductance peak if this splitting is larger than the Andreev level 
broadening (proportional to t r ) caused by the coupling of the ferromagnet to the reservoir through a potential barrier 
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of the transparency t r <C 1. In Fig. [5] this splitting is absent at tp = tt because in this case <5rjv = 0.05 < t r = 0.1. 
However, with a deviation of the superconductor phase difference tp from tt the splitting sharply increases, and a dip 
in the conductance peak appears as is seen in Fig[S]at ip = \.\tt and p = 1.2n. 

We note here that the differential conductance G as a function of eV is proportional to the the density of Andreev 
states in the ferromagnet permitting a direct spectroscopy of Andreev levels by conductance and current measure- 
ments. 



DENSITY OF STATES IN DIFFUSIVE S/F/S STRUCTURES WITH I < |A| 

As is shown in Section „ for ballistic S-F-S structures, the energy level £ — Iq < |A| is Aj_-fold degenerated and the 
density of states v{e) at e = Iq and ip = tt is proportional to N\ ~ S/\ 2 F 3> 1 (see Fig. 0} in the same manner as it 
take place in S-N-S (non- magnetic) structures at e = (see |E^,|Hf|). In diffusive non-magnetic S-N-S structures there 
is an energy gap in the density of states around e = which is maximal at p = and shrinking to zero together with 
a sharp increase of the density of states at e = as tp approaches tt . In order to see whether the above-mentioned 
degeneracy of the level e = Iq < |A survives in diffusive S/F/S structures, we consider the density of states in a 
diffusive S/F/S structure in the vicinity of tp = n(2k + 1), k = 0, ±1, ±2, ... assuming the motion of quasiparticles 
inside the ferromagnet to be semiclassical. It allows us to use Gutzwiller's approach [§4| which has the advantage of the 
clarity of physical presentation. On the other hand, the new class of two-dimensional (2D) magnetic semiconductors 
with large dielectric constants and small effective masses [§2 very well satisfies the condition of semiclassical motion 
a — r s pp/Ti ;§> 1 (r s is the screening length in the ferromagnet) with a = 10 -f- 10 2 . We show, that in much the same 
way as in the ballistic case (see Section B and Section^, in a diffusive S/F/S structure the energy level £ = Iq (Iq < |A|) 
is also Nq ~ S/ \ 2 F -fold degenerated and the hight of the peak in the density of states at e = Iq is proportional to 
Ao, if the phase difference between the superconductors is in the vicinity of an odd number of tt. Below we consider 
the case that the Andreev reflection of a quasiparticle at the F/S interfaces is accompanied by a normal specular 

(1 2) 

reflection in the same manner as for the ballistic case, that is the probability amplitudes \r N ' | <C 1. 

For calculations of the density of states ^sfs(s) at p = n(2k + 1), k = 0, ±1, ... we use the Gutzwiller semiclassical 
trace formula which shows that ^sfs(^) can be presented as a sum over all periodic classical trajectories |6^| : 



v S Fs{s)=Y. A ^y S3[e) ( 27 ) 

j 

(here j is the trajectory index, Sj — § pdq is the classical action integral (the integral is taken along the periodic 
j-trajectory), Aj is the trajectory amplitude which is equal to the Gaussian path integral around the stationary phase 
path (which is the classical trajectory)). 

If there are Andreev and normal reflections at the F/S boundaries, the semiclassical trajectories contributing to 
the density of states Eq. H27J1 are of the type shown in Fig. In this case the problem of finding the density of states 
is reduced to a quantum scattering problem for the configuration shown in Fig. 03. The points of reflections at the 
F/S interfaces are shown with black dots. Propagation between these points is coherent in both electron and hole 
channels, which is illustrated by dashed and solid lines of equal lengths. For a quasiparticle with a non-zero energy 
£ the phase gains along electron and hole trajectories do not completely cancel since the momenta are now different. 
The resulting decompensation effect is of order of ||e| — Iq\/Eth [SSj 

( 1 21 

If r N ' ^ 0, the motion of a quasiparticle along such a trajectory is of a quantum character despite the quasiclassical 
parameter is small because quasi-particle waves split at the scattering points. The semiclassical wave- functions in the 
adjacent sections of Fig. [7}3 are connected by the 2-channel scattering matrix [7| describing Andreev and normal 
reflections at the F/S boundaries. The problem of finding the spectrum and the wave- function in such a one- 
dimensional chain is reduced to solving a set of matching equations for amplitudes of electron off and hole <z„ 

(12) 

excitations in every section of propagation between scattering points. If r N ' <C 1, the main contribution to the 
density of states Eq.J23 at tp = 7r(2fc + 1), k = 0, ±1, ±2, ... comes from trajectories which do not include successive 
reflections at the same F/S boundary. Taking these observation into account the following set of algebraic matching 
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FIG. 7: A, Structure analyzed in the text with superconducting (S) and disordered ferromagnetic (F) elements labeled. The 
ferromagnetic element is coupled to two superconductors. Semiclassical trajectories of electron- and hole-like guasiparticle 
excitations are indicated by solid and dashed lines, respectively. Sections of the trajectories going between F/S boundaries 1 
and 2 are connected by Andreev and normal reflections at the points shown as black dots. B, Sequence of scattering events 
along the trajectory shown in A; L n is the length of the trajectory in section n. 



equations emerge: 

(e) _ (2) »W (e) (2) i9 W (h) 
a JV a -l + r A e a 

a W = _ r P)* e i*L1 a W +r P)* e _,*w a (/i) 



4 A) = -r^e^o'a^ + r^e^a^ (28) 



(e) (1)* (e) , (1)* (h) 

a_2 — ~ r A e a -2+ r N e '"-1 



The amount of the phase gain after propagation along the trajectories in section n is 

$ " = /l P " ' H)dl,n * $ " 0) ± T ™ (£ ~ aI ° )ln 

where 4>^°- ) = ppL n /h and r„ = L n /vF is the propagation time in section n. The quantity r^ L,2 ' ) and r^' 2 ^ are, 
respectively, the probability amplitudes for Andreev and normal reflections at F/S boundaries 1 and 2. Phases and 
amplitudes of electrons or holes along semiclassical paths arc defined in such a way that no phase has been gained at 
the beginning of a particular electron or hole section n. Hence the amplitude is al 6 '' 1 ' at the beginning of the section, 
and oil ^ exp ±$£ e '^ at its end. We note that one can show that the large phases <&i°' ) can be removed from the 
set of eqns (|28|) . This is a manifestation of the fact that the electron and hole phase gains compensate each other at 

£ = al . 

The set of equations Ea. (|28(l can be presented in a matrix form: 



|a>=C/|a> 



(29) 
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where components of the vector |a > are the coefficients ajf of Ea (|28() : the definition of the matrix U is obvious 
from Eq.(|2HJ: it is a unitary matrix written as follows: 

U mk (& e >») = V mk exp{i^ h) } (30) 

where the m-th row of the matrix V has only two non-zero elements which are elements of the matrices T^ 1,2 ) (see 

Eq.Q); *(^) = (...,<& { ^\<S>[ e ' h)M -' h) ,...). 

According to general principles of quantum mechanics |60| , the spectrum of the system is defined by zeros of the 
determinant 



Det 



I -U(*(e)) =0 (31) 



where / is the unit matrix. 

Summation in Eq. (|27|l is over different trajectories and one can think of it as averaging over various distributions 
of r n . Therefore, the density of states Ea.l(2*7)l for the S/F/S case under consideration can be re- written as follows: 



vsfs 



(e)^N J2(Wrand Om m (32) 



where N = S/X 2 F and << ... >> implies an averaging over r„; v° andom {e) is the density of Andreev states for a 
fixed spin projection a generated by a given semiclassical trajectory of the type shown in Fig. UJ3. The distribution 
of propagation times r„ depends on details of the disordered potential in the mesoscopic ferromagnetic region. These 
are not known, but it is natural to assume that the propagation times along different sections of the semiclassical 
trajectories (see Fig. |7JV) are uncorrected. 

In order to find the averaged density of states << v° andom {e) >> we use the presentation suggested by Slutskin 
|6ll | for the density of states of a one-dimensional chain of the type in Fig. 03. In the general case it can be written 
as follows: 



V( £ ) = ^J2 T ^ F ™ +F nn + 1 ) (33) 
n=l 

where matrix F satisfies the matrix equation F = FU + U (in our case U is defined by Ea. ^OjO and 

F nn = ^A„(l)exp{ilf} (34) 
l 

Here I = {h, ■■ -,In} (N is the number of sections in the chain), l n are integers, either positive or equal to zero, 
If = Ylk=i Ikfki f — {/i( £ )i •••) /iv( £ )}> /n( £ ) is the phase gain along section n of the chain, Fourier coefficients A n 

(1 2) 

depends on "hopping integrals" r N ' and are independent of t^, k — 1,2, ...N; summation is over integer l k , either 
positive or equal to zero. Therefore, the averaged density of states is 



1 N 

(Wandomie))) = E ^ + F « + W ^ 



n=l 



where F n = F nn . 

For the case under consideration fk — Tk(e — (tIq)/% and for the convenience sake we re-write Ea. H34|) as follows: 

F n = ]T A n (l) n exp {il k r k (e - al )/h} (36) 
l fe 

Since amplitudes A n do not depend on r k one only has to average the product T n Y[ k expil k T k (e — crIo)/h while 
calculating the average density of states Eq. . 

For the configuration of Fig. [7)3 the averaged density of states can be found exactly if one chooses a Lorentz form 
for the distribution function P(t), 
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Using the Lorenzian distribution Ea. (|37|l one finds the result 

-2- / d £l 

J— (ei - (e - ^o)) z + (7/^) 2 (e - ^o) 2 



t/c ^QexpiZ fe Tfe(e - al )/h \\ = / cfci- 2 £l — I r JJexpiZ/-T£i/7i j (38) 



Inserting Eq. Ij38(l in Eq. (|35J) one finds the averaged density of states 

4 + 4( 7 /f) 4 (e-aIo) 

where 



1 N 



n=l 



AT 

/» V 
and 

F n = £ 4.00 II ex P {^/M ( 41 ) 
l fe 

While writing Ea. (|39|l we took into account the fact that v p {e\) — v p (—ei). 

From Ea. (|40|l and Eq. (|41|l one sees v v (e) to be the density of states for the case that the system in Fig. [7)3 is 
non-magnetic (Jo = 0) and periodic with all r n = f (that is all L n = L where L =<< L n >>= t/vf ~ (vf/D)L 2 
where L is the distance between the S/F interfaces). 

For the case of a periodic chain, Ea. (|28() with all t„ = f can be easily solved, and for the phase difference ip in the 
vicinity of an odd number of tt one gets the dispersion low for the quasiparticle moving along the chain as follows: 



e± (*) = ±E Th V(M 2 + \r { N ] I 2 + l^ 2) I 2 ~ 2|r# | |r£> | cos fc (42) 

where iJy^ = hD/L 2 is the Thouless energy, Stp = (p — lit (I is an odd number, |<$^| <C 1), A: is a continuous quantum 
number - the " quasi- momentum" in the periodic one-dimensional chain. 
Using Eq. (|42fl one finds the density of states v v of the periodic chain to be 

1\f\ 

vJe) = 11 (43) 

wJ(e 2 - e 2 ■ ) (e 2 ^e 2 ) 

"V \ ^minJ \ max > 

where the minimal and the maximal energies of the energy band Eq. (|42|1 are 

(44) 



)i _ | r ( 2 )| N2 

\'N 



E Th ^j{5^) 2 + (\r^ 
-E Th \j{5 V ) 2 + (\r^\ + \rf\) 



2 



Inserting Eq . (|44|l in Eo. (|39|l and integrating the resulting expression, we find an exact formula for the averaged 
density of states of the random chain in Fig. [7J 

1/2 



2k\s-ctI \ \ ^{e A mm + 4(«|e - tr/pl) 4 ) (e 4 ^ + 4(«|e - a/ |) 4 ) + e 2 min e 2 max - 4(«|e - aJ |) 4 

random v / S rj / 4 



2 (4«„ + 4(«|e - ^o|) 4 ) (e 4 „a, + 4(«|e - ^ol) 4 ) j 

(45) 

where k = 7/1". 

Multiplying Ea. H45|) by N and summing over the spin projections (see Ea. (|32[l ) we get the average density of states 
of a diffusive S-F-S structure as follows: 



v 2\e -oIq\ f V(ef mra + 4(e - a/ |) 4 ) (g^ + 4(e - <t/ |) 4 ) + e 2 mm e 2 max - 4(e - a/p) 4 1 1/2 
S ^ SlJ °V 2(e 4 im + 4( £ - ( T/ |) 4 )( £ 4 „ a:c + 4( £ -a/o|) 4 ) f 



(46) 



(we have used 7 = f, see Ref. [56| ) 

Numerical results for the averaged density of states for a diffusive S/F/S structure based on Eq . 1^461) are presented 
in Fig. |SJ The sharp peak in the dependence of the averaged density of states on energy In Fig. |SJ3 corresponds to 
the iV -fold degeneracy of the energy level e = Iq at <p equal to odd numbers of 7r; the splitting of the peak in Fig. 
IHK is proportional to \\r$ \ — 
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A B 

FIG. 8: A, The averaged density of states for a diffusive S/F/S structure for To = Etk and Tn = 0.05, r^f 1 — 0.1, ip = 7r. B, 
The averaged density of states for a diffusive S/F/S structure for Iq = Etk and r£ = ri 2 'o.l; fo = 2/(jTETh) 



CONCLUSION 



In conclusion, for S/F/S structures with I < A (I is the ferromagnet exchange energy, A is the superconductor 
energy gap) we have demonstrated that an extremely high degeneracy of the energy level e = Iq proportional to 
No = S/\ 2 F at ip equal to odd numbers of ir (S is the cross-section area of the ferromagnet, Xp is the electron 
wave length) results in giant oscillations of the density of Andreev states and the conductance of the ferromagnet 
with a change of the superconductor phase difference <p. This phenomena is a convenient tool for the Andreev level 
spectroscopy (the differential conductance G is proportional to the density of Andreev states) and enables applications, 
e. g. as a double-gate ferromagnet transistor analogous to the one described in 50] . On the other hand, this effect 
permits to find the exchange energy, Iq, of the interaction between the electron spin and the spontaneous moment 
of the ferromagnetic conductor by an electric measurement of the differential conductance because a sharp and giant 
peak in the dependence of the conductance on the applied voltage V arises at V = 2Iq / e (the corresponding peak in 
the density of Andreev states takes place at energy e = Iq; see Fig0] Fig. Eland Fig. |SJ|. 
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